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Electronic  Factors  in  Thiophene  Adsorption  and  Hydrodesulfurization 
on  M0S2  Surfaces 

MARJANNE  C.  ZONNEVVLLE  and  ROALD  HOFFMANN 

Department  of  Chemistry  and  Materials  Science  Center,  Cornell  University,  Ithaca,  New 
York  14853,  U.S.A. 

Hydrodesulfurization,  the  removal  of  sulfur  in  the  form  of  H^S  from 
petroleum  compounds,  is  a  crucial  step  in  industrial  refinement  processes. 

In  spite  of  numerous  recent  experimental  and  theoretical  studies,  many 
aspects  of  the  complex  HDS  system  remain  obscure-* .  We  will  focus  on  the 
MoSg/thiophene  system  in  this  work. 

M0S2  based  systems  are  among  the  most  frequently  used  and  studied 
catalysts.  The  microscopic  natur/of  the  active  surface  is  not  well 
defined,  although  several  points  can  be  made.  Crystalites  of  the  material 
are  bound  to  an  AI2O3  support,  which  is  generally  considered  to  be  inert. 

M0S2  is  a  layer  compound;  the  crystals  cleave  readily  between  the  S-Mo-S 
sandwiches  to  expose  all-sulfur  faces.  However,  this  surface,  the  basal 
plane,  is  believed  to  be  inactive  in  HDS.  The  catalytic  surface  contains 
exposed  Mo  atoms.  A  good  candidate  is  the  edge  plane  perpendicular  to  the 
sulfur  layers.  Corner  sites  along  the  edge  and  defect  sites  arising  from 
"missing"  surface  sulfur  atoms,  in  particular,  have  been  implicated  as  the 
active  sites2.  These  Mo  atoms  are  coordinatively  unsaturated,  thus  it  is 
not  unreasonable  that  adsorption  may  occur  here,  u:-  .  ’ 

Thiophene  (SC4H4)  is  the  simplest,  and  most  reactive,  aromatic 
species  in  the  HDS  feedstock.  Neither  the  geometry  of  adsorption  nor  of  the 
active  site  is  known.  Coordination  of  the  planar  five-membered  ring 
species  directly  to  a  surface  Mo  is  generally  assumed  for  both.  In  organo- 
metallic  chemistry,  both  rj5-bound  (coordinated  to  all  five  atoms  of  the 
ring)  and  tj1  -sulfur  bound  (coordinated  perpendicularly  to  the  sulfur) 
species  exist^.  However,  on  the  M0S2  surface,  r|2  coordination  via  the 
p-carbons,  which  are  furthest  from  the  sulfur,  or  via  the  S-C  bond  cannot 
be  excluded.  The  severe  experimental  difficulties  associated  with  a  direct 
analysis  of  the  reaction  mechanism  on  the  catalytically  important  surfaces 
have  been  prohibitive.  However,  model  studies  have  been  performed  on 
various  transition  metal  surfaces.  On  Mo(100)  and  Ni(100),  evidence  for  a 
metallocycle  formation  has  been  found4,  whereas  on  Mo(1 10)  the  proposed 
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intermediate  is  an  ri1-coodinated  surface  butyl  thiolate^. 

As  the  organic  products  are  mainly  butanes  and  butenes,  the  S-C  bonds 
must  be  severed  in  any  case.  Bonding  modes  which  induce  partial  occu¬ 
pation  of  the  thiophene  LUMO,  3bi ,  will  weaken  the  bond.  It  is  a  n*  orbital 
antibonding  between  S  and  C,  and  localized  mainly  on  the  a-carbon  (adja¬ 
cent  to  the  sulfur).  The  extent  of  3bi  occupation,  as  measured  by  the  inte¬ 
grated  projected  density  of  states,  is  one  criteria  we  can  use  to  identify 
possible  active  sites.  The  orbital  lies  4.4eV  above  the  HOMO,  thus  sites 
which  produce  substantial  occupation  are  also  likely  to  be  high  energy 
sites. 


To  simplify  our  examination  of  the  HDS  active  sites,  our  extended 
Huckel  calculations  were  performed  on  a  three  layer,  one-dimensional 
infinite  chain  M0S2  sandwich  such  as  1  (shown  in  with  one  particular 
thiophene  coordination  mode).  The  geometric  specifications  are  taken 
from  the  real  crystal  strucure.  Interadsorbate  interactions  were  mini¬ 
mized  by  using  a  unit  cell  containing  six  atoms,  four  sulfur  and  two 
molybdenum,  per  layer.  One  such  cell  is  shown;  x  is  the  direction  of 
propegation.  Atomic  charges  agree  to  within  0.005e-  and  overlap  pop¬ 
ulations  to  0.005  with  those  calculated  for  a  two-dimensional  M0S2  slab 
(propegated  in  both  x  and  y  directions)  with  the  bulk  geometric  parameters. 
The  one  dimensional  model  is  appropriate  for  this  problem  as  "inter¬ 
sandwich"  interactions  are  negligible  (at  the  extended  Huckel  level)  due  to 
the  3.1 7  A  spacing  between  the  sulfur  layers.  Similar  agreement  in  the 
calculated  parameters  is  found  between  one-  and  two-dimensional  slabs 
with  thiophene  adsorbed  at  several  test  sites.  The  atomic  parameters  used 
can  be  found  in  reference  6. 

Our  model  exposes  two  types  of  surface  Mo  atoms.  M03  protrudes  from 
the  surface  and,  with  only  four  neighboring  sulfur  atoms,  is  coordinatively 
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unsaturated.  Mog  is  recessed  into  the  face  and  has  the  full  bulk  compli 
-ment  of  six  neighboring  sulfurs.  Both  are  electron-deficient,  even  with 
respect  to  the  bulk  (Mo:  +.66,  S:  -.33);  the  net  charges  are  +0.80  for  M03  and 
+1.13  for  Mog.  Seven  geometries  were  considered  for  thiophene  adsorption 
onto  a  pristine  face  (no  defect  sites).  For  the  first  set  the  thiophene 
sulfur,  S{,  is  bound  on  top  of  M03  such  that:  1)  the  ring  is  perpendicular  to 
the  surface  in  the  xz  plane,  2)  the  ring  is  at  45°  to  the  surface,  and  3)  the 
ring  is  parallel  to  the  surface.  The  same  ring  orientations  are  selected  for 
4,  5  and  6  such  that  is  bridging  M03  and  Mog.  The  sixth  is  an  ti5  species 
and  is  shown  in  the  schematic.  The  coordination  of  7  is  as  for  1,  but  to  the 
receding  Mog.  The  St-Mo  distance  is  1 .90A  for  all  modes,  which  is  the 
minimum  found  in  Mo  organometallic  species^.  Similar  surface-adsorbate 
calculations  indicate 

that  the  choice  of  chemisorptive  bond  length  will  alter  the  magnitude, 
rather  than  the  direction  and  type,  of  interaction8. 

On  the  basis  of  both  the  S-C  overlap  population  (o.p.)  and  the  S-C 
antibonding  3b-|  occupation,  the  bridging  sites  are  substantially  more 
effective  at  weakening  the  S-C  bond  than  the  on-top  modes  in  any  given 
ring  orientation.  Coordination  to  the  recessed  Mo  has  the  smallest  effect 
on  this  bond.  The  i\ 5  species,  6,  experiences  the  most  severe  S-C  bond 
weakening;  the  o.p.  drops  from  the  molecular  value  of  0.887  to  0.778.  The 
thiophene  3b-|  occupation  of  0.71 9e-  is  the  largest  calculated.  By 
matching  the  peaks  in  the  3b-|  DOS  to  those  of  the  surface  Mo  levels,  we 
find  the  primary  source  of  backdonation  to  be  the  Mog  x2-y2  levels.  The  xy 
states  are  involved  in  the  M0S2  surface  bond  formation;  ~40%  are  found  in 
the  sulfur  p  block.  In  a  true  Cgv  fragment,  CpM  (Cp=cyclopentadienyl, 
C5H5),  the  xy  and  x2-y2  orbitals  are  degenerate  and  both  interact  with  the 
degenerate  Cp  k*  orbitals8,  which  become  the  3b-|  and  2a-|  of  thiophene.  If 
sulfur  atoms  are  sequentially  removed  from  the  surface  to  model  defect 
sites,  the  xy  levels  shift  out  of  the  sulfur  p  block  and  into  major  peaks  of 
the  thiophene  n*  states. 

The  3b-|  occupation  increases  for  the  bridging  set  in  the  order  perp¬ 
endicular,  45°  and  parallel,  but  much  less  so  for  the  on-top  set.  As 
described  earlier,  the  binding  energies  become  less  favorable  as  well.  One 
must  keep  in  mind,  however,  that  our  calculations  produce  reasonable 
relative  energies,  but  unreliable  absolute  energies. 

Chemisorptive  bonds  typically  have  two  components:  adsorbate  to 
surface  donation  and  surface  to  adsorbate  backdonation.  The  two  sulfur 
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lone  pair  orbitals,  2b-|  perpendicular  to  the  ring  and  9a-|  in  the  ring  plane, 
figure  prominently  in  the  former.  Thiophene  is  a  much  better  donor  than 
acceptor  adsorbate,  as  the  LUMO  is  not  only  localized  at  the  a-carbon  and 
very  high  in  energy,  but  the  sulfur  lone  pair  orbitals  are  only  ~2eV  below 
the  fermi  level.  The  projected  DOS's  indicate  that  for  each  of  the  geo¬ 
metries,  the  more  important  of  the  two  is  the  one  best  able  to  interact 
with  the  Mo  yz  states.  For  the  on-top  geometries,  the  interaction  is 
jr-bonding,  involving  2b-|  for  the  perpendicular  site  1  and  9a-j  for  the 
parallel  3.  The  lone  pair  oriented  radially,  rather  than  tangentially,  with 
respect  to  the  Mo-St-Mo  angle  dominates  the  bridging  modes  ( 9a-)  for  the 
perpendicular  4  and  2b-|  for  the  parallel  6).  A  comparison  of  the  dis¬ 
persion  and  charge  transfer  out  of  the  lone  pairs  shows  that  bridging 
thiophene  is  a  superior  electron  donor  to  on-top  thiophene.  The  direction  of 
net  charge  flow  is  exclusively  from  the  adsorbate  to  the  substrate,  and  into 
the  uncoordinated  surface  or  bulk  Mo  atoms.  Coordinated  Mo  atoms  tend  to 
be  further  oxidized.  The  most  profound  effect  is  the  +1 .27  charge  at  the  r|5 
Mo  relative  to  the  bare  surface. 

The  sequential  removal  of  three  of  the  four  sulfur  neighboring  a 
coordinated  Mo  atom  brings  out  a  problem  inherent  in  non-infinite  models, 
that  being,  what  charge  to  assign  to  the  non-stoichr  metric  unit  cell.  The 
problem  will  be  addressed  in  more  detail  in  a  subsequent  work.  Regardless, 
the  effect  on  the  dispersion  and  distribution  of  the  thiophene  orbitals  is 
significant  in  most  coordination  geometries.  We  discuss  only  the  r|5  case 
here.  As  previously  mentioned,  the  Mo  xy  becomes  increasingly  involved  in 
backbonding.  The  important  thiophene  3b-j  and  2b-|  become  more  compact 
as  they  interact  with  less  difuse  Mo  orbitals  (fewer  M0S2  surface  bonds). 
Assuming  a  constant  fermi  energy  to  circumvent  the  aforementioned  charge 
problem  leads  to  a  slight  increase  in  the  3b-(  occupation  (-5%  total)  and 
decrease  in  the  S-C  o.p.  (-10%  total). 

Our  extended  Huckel  calculations  implicate  r^-bound  thiophene  as  a 
possible  active  species  for  HDS.  Defect  sites  may  be  more  reactive. 
Certainly  more  geometries  must  be  explored  before  a  definitive  active  site 
is  discovered. 

Acknowledgement.  We  wish  to  thank  Suzanne  Harris  for  helpful 
disscusions,  and  the  Office  of  Naval  Research  for  their  support  of  this 
work. 

References 

la.  Zdrazil.  M..  AddI.  Catal.  4.  107  M9821. 


5 

b.  Massoth,  F.E.,  Advan.  Catal.  27,  265  (1 978). 

c.  Harris,  S.,  Chianelli,  R.R.,  J.  Catalysis  86,  400  (1 984). 

2a.  Roxolo,  C.B.,  Daage,  M.,  Ruppert,  A.F.,  Chianelli,  R.R.,  J.  Catalysis  100, 

176  (1986). 

b.  Joffre,  J.,  Geneste,  P.,  Lerner,  D.A.,  J,  Catalysis  92,  56  (1 985). 

3a.  Lesch,  D.A.,  Angelici,  R.J.,  J.  Amer.  Chem.  Soc.  106,  2901  (1984). 

b.  Kuehn,  C.A.,  Taube,  H.,  J.  Amer.  Chem.  Soc.  98,  689  (1976). 

4a.  Zaera,  F.,  Kollin,  E.B.,  Gland,  J.L.,  to  be  published, 
b.  Zaera,  F.,  Kollin,  E.B.,  Gland,  J.L.,  to  be  published. 

5.  Roberts,  J.T.,  Friend,  C.M.,  to  be  published. 

6.  Summerville,  R.H.,  Hoffmann,  R.,  J.  Amer.  Chem.  Soc.  98.  7240  (1976). 

7.  Bart,  J.C.J,  Ragaini,  V.,  in  "Proceedings  of  the  Climax  Fourth 
International  Conference  on  the  Chemistry  of  Molybdenum"  (H.F.  Gray, 

P.C.H.  Mitchell,  Eds.)  pi  9.  Climax  Molybdenum  Co.,  Ann  Arbor,  1979. 

8.  Zonnevylle,  M.C.,  Hoffmann,  R.,  to  be  published. 

9.  Albright,  T.A.,  Burdett,  J.K.,  Whangbo,  M.H.,  "Orbital  Interactions  in 
Chemistry".  John  Wiley,  1985,  New  York,  p.388. 


'•..n  i  or 


;  r4  *r 

l  -  - 


.U: :  c 


£ 

□ 


ju.iuil'ieritioa. 


I  TV 


'  f; T. v iX’.t.  i on/ 


JY VO \  1 1 . L- i  1  i t y  c o d a 3_ 
_ Avail  and/or 


X 


